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ABSTRACT 

Using the radio observations by FIRST and NVSS, we build a sample of 151 
radio variable quasars selected from the Sloan Digital Sky Survey Data Release 
3 (SDSS DR3). Six (probably another two) among them are classified as broad 
absorption line (BAL) quasars, with radio flux variations of a few 10 percent 
within 1.5-5 years. Such large amplitudes of the variations imply brightness tem- 
peratures much higher than the inverse Compton limits (10^^ K) in all the BAL 
quasars, suggesting the presence of relativistic jets beaming toward the observer. 
The angle between the outflow and the jet is constrained to be less than ~ 20°. 
Such BAL quasars with polar outflows are beyond the simple unification models 
of BAL quasars and non-BAL quasars, which hypothesize that BAL quasars are 
normal quasars seen nearly edge-on. 

Subject headings: quasars: general — quasars: absorption lines — radiation: 
radio continuum 



1. Introduction 

The broad absorption lines (BALs), which appear in the UV spectra of ~ 15% quasars 
(e.g., Tolea et al. 2002; Reichard et al. 2003a; Hewett & Foltz 2003), are characterized by 
prominent broad and blueshifted absorption troughs due to ions of a wide range species from 
Fell, Mgll to OVL It is now commonly accepted that BAL region is present in all quasars 
but with covering factor much less than unit. The dichotomy of BAL quasars and non-BAL 
quasars is interpreted as an orientation effect. For instance, Murray et al. (1995) suggested 
that BALs present themselves in a minority of the quasar population when the line of sight 
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passes through the accretion disk wind nearly at the equatorial plane. A detail study shows 
that such a scenario is also consistent with the continuum polarization and X-ray absorption 
in BAL quasars (Wang, Wang & Wang 2005). Elvis (2000) appealed to a funnel-shaped thin 
shell outflow that arises from the accretion disk to explain various observational properties 
of quasars, including BALs which are interpreted as "normal" quasars when the wind is 
viewed end-on. Both orientation models require a rather large incline angle of i ~ 60° for 
BAL quasars. 

The determination of the inclination of the accretion disk in BAL quasars is vital to 
understand the geometry of the BAL outflow. The axis dcflncd by rclativistic radio jets, 
which is likely aligned with that of the accretion disk (Wills et al. 1999), can be used to 
infer the inclination angle of the accretion disk. Becker et al. (2000) found that the 29 
BAL quasars discovered in the FIRST Bright Quasar Survey (FBQS) exhibit compact radio 
morphologies (~ 80% unresolved at 0".2 resolution) and show wide scatter in the radio 
spectral index (a ~ —0.7 — 1.2 and a < 0.5 for ~ 1/3 sources, S^, oc u~"). According to 
the unification model of radio-loud AGNs (cf. Urry & Padovani 1995), core-dominated flat- 
spectrum radio sources are those viewed close to the radio jet axis, while lobe-dominated 
steep-spectrum radio sources appear to be present at larger viewing angles. The radio 
morphology and spectra of FBQS BAL quasars indicate that the orientation of BAL quasars 
can be both face-on and edge-on, contrary to the simple unification model introduced above. 
However, their radio spectral indexes, derived from non-simultaneous observations, might be 
biased by radio variations. In addition, most FIRST detected BAL QSOs are only radio- 
intermediate, for which it is unknown if the unification model based on radio-loud AGN (i.e., 
core-dominated fiat-spectrum radio sources are face-on) still applies. Also the size of radio- 
intermediate sources might be much smaller, such as the radio intermediate quasar 111 ZW 
2 (Brunthaler et al. 2005), thus observations with much higher spatial resolution is required 
to confirm its compactness. Jiang & Wang (2003) acquired high resolution VLBl images of 
three BAL quasars at 1.6 GHz. They found one source is resolved into asymmetric two-sided 
structure and a bright central component at ~ 20 mas resolution. This morphology mimics 
a Compact Steep Spectral (CSS) source, but its size is much smaller than the typical value 
for a CSS source. The other two sources remain unresolved at this resolution, indicating 
that they are viewed face-on. 

The radio fiux variability, which is sometimes a better indicator of jet orientation than 
the radio morphology (at a limited resolution and sensitivity) and the spectral slope (often 
based on non-simultaneous observations), has not hitherto been adequately explored. Becker 
et al. (2000) commented that 5 BAL quasars in their sample appear to be variable at 1.4 
GHz, but no further concern was given. In this paper we use two-epoch observations at 
the same frequency of 1.4 GHz by FIRST (the Faint Images of the Radio Sky at Twenty 
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centimeters survey, Becker et al. 1995) and NVSS (the NRAO VLA Sky Survey, Condon et al. 
1998) to study the radio variabihty of the BAL quasars observed by SDSS (the Sloan Digital 
Sky Survey, York et al. 2000). We identify 6 (probably another two) BAL quasars with 
radio variation with a few 10 percent. Calculations based on the radio variations imply that 
these sources should be viewed face-on with inclination angles less than 20°. This confirms 
the existence of polar BAL outflows, opposite to the unification model for BAL quasars. 
Throughout the paper, we assume a cosmology with Hq~ 70 km Mpc^^, 0,m = 0.3, and 
Qa = 0.7. 



2. Data and Analysis 

2.1. The SDSS Quasar Catalog and the Radio Data 

Our jumping-off point is the SDSS Quasar Catalog (the 3rd edition, Schneider et al. 
2005, hereafter S05C for the catalog). Its sky coverage is ~ 4, 188 deg'^, which accounts for 
more than 2/5 of ~ 10^ deg'^ planed by SDSS. S05C consists of the 46,420 quasars with 
iband absolute magnitudes Mi < —22.0 and positional uncertainties < 0".2. In particular, 
optically unresolved objects brighter than i — 19"*. 1 with FIRST counterpart within 2" are 
observed, which we are more interested in. 

Using the NARO Very Large Array in its B configuration, the FIRST survey began in 
1993 was designed to explore the faint radio sky down to a limit of Si,4^ghz ~ 1 rnJy with 
a resolution of ~ 5" FWHM. Its sky coverage is mostly superposed with that of SDSS. The 
FIRST radio catalog was presented in White et al. (1997) with a positional accuracy of < 1" 
(with 90% confidence level). By matching the SC05 and the FIRST catalog, we find 3,757 
quasars in S05C with FIRST counterparts within 2". Considering the average source density 
of ~ 90 deg'"^ in the FIRST survey, we expect only ~ 0.1% of these matches to be spurious. 
We note that the above cutoff is biased against quasars with extended radio morphologies, 
e.g., lobe-dominated quasars. However, missing this kind of radio variable sources based on 
the data presently available (see §2.2) and this does not influence our main results. After 
all, only ~ 8% lobe-dominated radio sources will be lost (Ivezic et al. 2002; Lu et al. in 
preparation), and such sources only accounts for ~ 3.1% in a well-defined variable sample of 
de Vries et al. (2004). 

NVSS, using the VLA in its more compact D and DnC configurations, was carried out 
between 1993—1996 at the same frequency as FIRST. NVSS covers all the sky area of FIRST, 
but with a higher survey limit of Si^ghz ~ 2.5 mJy and a lower resolution of FWHM ~ 45". 
The positional uncertainties are estimated to be varying from < l" for bright sources with 
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flux density > 15 mJy to 7" at the survey limit. With a typical background noise of 0.45 
mJy (about 3 times higher than that of FIRST), we expect that the NVSS survey be able to 
detect all of the FIRST sources with flux density Siaghz ^ 5 mJy, provided that the radio 
sources do not variate. 



2.2. The Selection of Radio Variable Quasars 

We first select quasars in S05C with redshift > 0.5 so that the Mgll absorption trough 
(or other bluer troughs), if presents, falls within the wavelength coverage of the SDSS spec- 
trograph (3,800 — 9,200 A). Out of these quasars, 1,877 have FIRST counterparts within 
2" and peak flux density Spp > 5 mJy as measured by FIRST. Then the FIRST images 
of all these 1,877 objects are visually inspected and classifled into three categories: 1) com- 
pact radio sources, 2) marginally resolved radio sources, and 3) radio sources with complex 
morphology. The first category includes 1,482 radio point sources unresolved at the FIRST 
resolution (79.0%). The second category includes 200 radio sources (10.7%), which often 
show core plus elongated (jet-like) structure, or a core embedded in weak diffuse emission. 
The third category includes 193 radio sources (10.7%), out of which 168 sources exhibit 
Fanaroff-Riley II morphology (FR-II, Fanaroff & Riley 1974). 

We search for NVSS counterparts for the quasars in all of the three categories within 
a 21" matching radius (3 a of the NVSS positional error at the survey hmit). NVSS coun- 
terparts are found for 1,838 of the 1,877 quasars within a 21" with a false rate of less than 
1%. Two possibilities may lead unfindable of the 39 quasars as NVSS counterparts within 
such a matching radius: 1) the flux falls below the NVSS limit resulting from variability; 2) 
the apparent centroid of the source is largely shifted due to contamination either by bright 
lobe(s) or by nearby unrelated bright source(s). We compare the NVSS and FIRST images 
to distinguish between the two and find that all of the 39 cases are due to confusion effects. 
In fact, 21 of the 39 quasars show FR-II morphology. These 39 quasars are excluded in 
our further analysis. This result indicates that care must be taken to compare the fluxes 
obtained during the two surveys. 

As noted by de Vries et al. (2004), the flux densities between FIRST and NVSS are not 
directly comparable since the radio sources that are not resolved by NVSS may be resolved 
by FIRST. In this case, the peak flux is smaller than the integrated flux density. We deflne 
the variabihty ratio (VR) for each quasar as 

VR = Sfp/Sni, (1) 
where Spp and S^j denote the peak flux density measured by FIRST and the integrated 
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flux density measured by NVSS. Our variability ratio is conservative for selecting variable 
sources which are brighter in the FIRST images. We estimate the signiflcance of radio flux 
variability as 

_ Spp — Sni /„s 

Cvar — / 2 ,2 
V ^FP + '^NI 

where ai\fj is the NVSS integrated flux uncertainty, and app the FIRST peak flux uncer- 
tainty. Objects with VR > 1 and ayar > 3 are taken as candidates radio variable quasars. 
This yields 154 candidate variable sources. We plot in Figure 1 (left panel) variabihty ratio 
VR against the peak flux density measured by FIRST Spp- The apparent dependence of 
VR on Spp is due to the combination of two facts, the dependence of measurement error on 
flux density and the confusion effects in the NVSS. The latter complication also induces the 
obvious asymmetric distribute of these sources around VR = 1. In fact, all but three (SDSS 
J160953.42+433411.5, SDSS J101754.85+470529.3, and SDSS J164928.87+304652.4) of ra- 
dio sources that are well resolved by FIRST have VR < 1. A more symmetric distribution 
can be found if only point radio sources and marginally resolved radio sources (symbols in 
blue and red color) are considered. Isolated compact sources located well below VR — 1 
are likely radio variable quasars. But new observations with higher resolution are needed 
to conflrm this. Out of the three quasars with complex radio morphology and VR > 1 (all 
are FR-II quasars), SDSS J160953.42+433411.5 has VR = 2.26 and a^ar = 20.90, fulflUing 
our selection criteria. After careful examination of its FIRST and NVSS image, we flnd 
confusion effects in NVSS are serious in this source and exclude it from our sample. Out 
of rest 153 candidates, 151 are point radio sources and 2 are marginally resolved. Their 
NVSS and FIRST images are all visually examined for possible contamination by unre- 
lated nearby bright sources, and one of them, SDSS J111344. 84-004411. 6, is removed for 
this reason. In addition, SDSS J094420. 44-1-613550.1 is eliminated from the sample because 
the NVSS pipehne gives a wrong flux density for this object. At last, we end up with a 
sample of 151 candidate radio variable quasars^. Using the FIRST integrated flux density 
and the SDSS PSF magnitudes, we calculate the radio-loudness (RL) of these quasars, de- 
fined as the k-corrected ratio of the 5 GHz radio fiux to the near ultraviolet flux at 2 500 A, 
RL — Si,^5GHz/S 25QQ^- ^ power law slope of a = {S^ a z/^") is assumed for radio emission 
since their radio spectra are likely flat (Wang et al. 2005), and the SDSS colors are used for 
the optical-UV emission. The radio properties of a sub-sample of the BAL quasars (see §2.3 
for detailed description) selected from these objects are listed in Table 1. In Figure 1 (right 
panel), we plot the radio-loudness against radio variability ratio for the 151 quasars. We see 
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that most radio sources with large variability amplitude {VR > 1.5) have the radio- loudness 
10 < RL < 250. 



2.3. The Optical Spectral Analysis and the Sample of Radio Vciriable BAL 

Quasars 

Eight BAL candidates were identified by visually inspecting the SDSS spectra of the 151 
radio variable quasars. These candidate spectra were corrected for Galactic extinction using 
the extinction curve of Schlegel et al. (1998) and brought to their rest frame before further 
analysis. We use the "Balnicity" Index (Bl) defined by Weymann et al. (1991) and Reichard 
(2003b) to quantitatively classify the absorption troughs. Following the procedures described 
by Reichard et al. (2003b), we calculated the Bis for the broad absorption troughs by 
comparing the observed spectra of our BAL candidates with the composite quasar spectrum 
created by Vanden Berk et al. (2001) from the SDSS Early Data Release (EDR). In brief, 
the EDR composite spectrum was reddened using the Pei (1992) SMC extinction curve to 
match each observed spectrum of our BAL candidates. The fit is done through minimization 
of with proper weights given to the wavelength regions that are obviously affected by 
prominent emission lines and absorption troughs. The results are displayed in Figure 2. 
Then we calculated the balnicity index of high ionization line BI(CIV) using the definition 
of Weymann et al. (1991), and the balnicity index of low ionization line BI(Mgll) and 
BI(AIIII) using the definition of Reichard (2003b): 

/•25,000 F''''^(lA 

BI = / dv[l - ,—^]C{v) (3) 



3,000 ' 0.9Ff^^{vy 

where F°'''{v) and Ff'*{v) are the observed and fitted fiuxes as a function of velocity in 
km s~^ from the systemic redshift within the range of each absorption trough, and 

^ _ I 1.0, if [1 — o^Qpfll ] > over a continous interval of > W km s~^, 
1 0, otherwise. 

The integral in Equation 3 begins at v — 3, 000 km for CIV and at = km for 
Mgll and AlIII. The given continuous interval in Equation 4 is 1^ = 2, 000 km s"^ for CIV 
and W — 1, 000 km for Mgll and AlIII. The results are listed in Table 1 and described 
individually below: 



SDSS J075310.42+210244.3 Apart from the deep high-ionization BAL troughs of CIV, 
SilV, NV with V ~ 0-13,500 km s~^, the low-ionization AlIII trough, covering a similar 
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velocity range, is also obvious. At the red end of the spectrum, the Mgll BAL trough 
is apparently present. With 5J(CIV)=3,633 km s"^ and S7(A1III)=1,420 km s~^, we 
are certainly observing a LoBAL quasar in this object. 

SDSS J081 102.91+500724.4 Based on the Balnicity Index of BI{CIV) = 617 km s'^ 
and the velocity range of ~ 6, 700 — 11, 600 km this object can be safely classified 
as a HiBAL quasar. Much shallower low-ionization BAL troughs of AlIII and Mgll 
with similar velocity range may also be present. However, high S/N spectrum is needed 
to confirm this. 

SDSS J082817.25+371853.7 Though its SDSS spectrum is rather noisy the low-ionization 
BAL troughs of AlIII and Mgll are securely identified. Classification of this object as 
LoBAL quasar should be safe based on our conservatively calculated BI of 1, 890 km 
for Aim and 626 km s'^ for MgIL 

SDSS J090552.40+025931.4 A BAL trough detached ~ 20, 000 km s'^ from the CIV 
peak is recognizable. However, the fit with the composite quasar spectrum yields 
BI{CIV) — km s~^. A marginal value of BI(CIV)=228 km s~^ is obtained if using 
a power law to fit the continuum. We tentatively classified this object as a candidate 
HiBAL quasar. 

SDSS J140126.15+520834.6 This object shows rather shallow BAL trough bluewards of 
CIV emission line. We label it as a candidate. 

SDSS J145926.33+493136.8 This object show two sets of BAL troughs, SilV. CIV. and 
Aim, around velocities v ~ 4,000 km s^^ and v ~ 15,000 km s^^. Wc classified it 
as a LoBAL solely by the presence of AlIII BAL troughs, since Mgll is redshifted 
out of the wavelength coverage of the SDSS spectrograph. Broad Lya emission line 
is completely absorbed and only narrow Lya emission line appears in the spectrum. 
Broad NV emission line is extremely strong perhaps due to scattering the Lya emission 
line (Wang, Wang, & Wang, in preparation). 

SDSS J153703.94+533219.9 We classified this object as a HiBAL quasar based on the 
presence of detached CIV BAL trough with BI{CIV) = 2, 060 km s'^. AlIII falls in 
the wavelength range of bad pixel and Mgll is redshifted out of spectroscopic coverage. 

SDSS J210757.67-062010.6 This object belong to the rare class of FeLoBAL quasars 
(e.g., Becker et al. 1997) characterized by the metastable Fell BAL troughs centered 
at 2,575 A and Mgll BAL troughs. Many Fell absorption features are also detected 
redward of Mgll, as well as neutral hehum absorption triplet, ife/AA2946, 3189, 3890, 
which can be used as powerful diagnostics of HI column density (Arav et al. 2001 and 
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references therein). Only a few quasars have been found to show Hel absorption hues 
(e.g., Anderson 1974; Arav et al. 2001; Hall et al. 2002). 



3. Discussion and Conclusion 

The origin of radio variability can be either extrinsic or intrinsic. The most familiar 
mechanism of extrinsic variability of radio sources is refractive InterStellar Scintillation (ISS, 
e.g., Blandford, Narayan, & Romani 1986). However, the variability amplitudes induced by 

ISS is seldom larger than ~ 2%. Considering the fact that all the 6 BAL quasars and the two 
BAL candidates we selected show radio variation with amplitude > 15% (Figure 1 and Table 
1), it is reasonable to believe that the radio variabilities of these BAL quasars are intrinsic 
to the radio sources. Marscher & Gear (1985) suggested that shocks propagating along the 
radio jets can induce flux variability. Amplitude of variability can be largely amplified by 
relativistic beaming effect if the radio jets are viewed nearly face-on. 

A lower limit of the brightness temperature can be inferred as follows (Krolik 1999), 

AP 

where AP^, is the variable part of the radio power computed from the difference between 
the FIRST and NVSS fluxes, At the time interval in the source rest frame between two 
observations, and ks the Boltzmann constant. We present Tg for the 6 BAL quasars and 
the 2 candidates in Table 1. We flnd that the brightness temperatures of all 8 radio sources 
are much larger than the inverse Compton limit of 10^^ K (Kellermann & Pauliny-Toth 1969). 
Such extremely large brightness temperatures strongly suggest the presence of relativistic 
jet beaming toward the observer. If the intrinsic brightness temperature of the radio sources 
is less than the inverse Compton limit, we can set a lower limit on their Doppler factor, 
— ( ioi2 Keivin V^^, hcucc au upper limit of the inclination angle, 

ei = arccos{[l-{'y6i)-']f3-'}, (6) 

where 7 = (1 — /3^)~^/^ is the Lorentz factor of the jets. We flnd that all these radio variable 
BAL quasars/candidates must be viewed within 6 < 10°, except SDSS J210757.67— 062010.6, 
for which 6 < 20°. If the equipartition value of ~ 5 x 10^'^ K (e.g., Readhead 1994; 
Lahteenmaki et al. 1999) instead of the inverse Compton value of 10^^ K is adopted as 
the maximum intrinsic brightness temperature, the inclination angle of our sample of BAL 
quasars should all be less than ~ 7°. Therefore polar BAL outflows must be present in 
these radio variable quasars, contrary to the simple unification models of BAL and non-BAL 
quasars, which hypothesize BAL quasars are normal quasars seen nearly edge-on. 
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Similar to most FIRST detected SDSS BAL quasars (Menou et al. 2001), All but one of 

our radio variable BAL quasars are all radio-intermediate with RL < 250. With RL — 923, 
SDSS J0828 17.25+371853. 7 is the only exceptionally very radio-loud BAL quasar in our 
sample. We note, however, its spectrum is significantly reddened with E{B — \/) ~ I*". 
Correcting for this intrinsic extinction would place the quasar in the radio-intermediate 
range with RL ~ 10. Since the radio emission is likely boosted greatly by the relativistic 
motion of the jet as indicated by its large brightness temperature, most of these sources 
would be intrinsically radio-weak or at most radio-intermediate. The properties of jets in 
radio-intermediate quasars are not well understood, but at least in one of such object III ZW 
2, super-luminal motion of knots at parsec scale has been detected (Brunthaler et al. 2005). 
Recurrent radio fiares was observed in the same object. Miller et al (1993) speculated that 
radio-intermediate quasars are actually relativistically boosted radio-quiet quasars based on 
the correlation between [0111] luminosity and radio power. It is possible that the BAL gas 
in these objects associated with the expansion of radio plasma. The origin of BAL in these 
objects may be different from the majority of BAL QSOs, in which a disk-wind is responsible 
for the absorption. 

Reichard et al. (2003a) found an uncorrected fraction of 14. Oil. 0% of BAL quasars in 
SDSS EDR within the redshift ranges 1.7 < 2; < 3.4, whereas Menou et al. (2001) found this 
fraction falls to 3. 3± 1.1% amongst those quasars with FIRST detection. Among the 82 radio 
variable QSOs within the same redshift range in our sample, we identify 4-6 BAL quasars. 
The overall fraction (~ 4%) of BAL quasars in our radio variable sample is similar to the 
fraction of the SDSS quasars detected by FIRST, but much lower (significant at 1-2(7 level) 
than the fraction of the whole SDSS quasars, the majority of which are selected according to 
optical color. As a comparison, out of ~ 600 quasars with z > 0.5 in the SDSS DR3 that are 
resolved by FIRST image, we only identified four BAL quasars. Their SDSS spectra with 
our model fits are displayed in Figure 3, and their FIRST images as contour maps. Three 
show FR-II morphology, out of which one object, SDSS J114111. 62-014306. 7 (also known 
as LBQS 1138-0126), was previously known as such (Brotherton et al. 2002). Apart from 
LBQS 1138-0126, the only other known FR-II BAL quasar is FIRST J101614.3+520916, 
which was also discovered from the FIRST survey (Gregg et al. 2000). At the sensitivity of 
the FIRST image, the majority of the resolved high redshift quasars are classical radio loud 
sources. The occurrence of BALs in such radio powerful quasars is extremely small (0.7%). 
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Table 1. Properties of the Radio Variable BAL Quasars 
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"Absorption trough falls out of the wavelength coverage of the SDSS spectrograph. 
^BI = 228 km if calculated using a power law fit to the continuum. 

Note. — Col.3: FIRST Peak fiux density, Col.4: FIRST Integrated flux density, Col.5: FIRST observation time, Col.6: NVSS Peak flux density, 
Col. 7: NVSS Integrated flux density. Col. 8: NVSS observation time. Col. 9: Significance of radio variability as defined in §2.2, Col. 10: lower limits 

of brightness temperature estimated as in §3, Col. 11: Radio-Loudncss calculated in §2.2, Col. 12 to Col. 14: Balnicity Index values of CIV, AlIII 
and Mgll, respectively. Col 15: our final classifications, question marked ones denote candidates. 
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Fig. 1. — The left panel shows the distribution of the 1,838 radio quasars on the flux density 
variation VR versus flux density Spp plan. VR is defined as as the ratio of peak flux density 
measured by FIRST (Spp) to integrated flux density measured by NVSS. Objects in the 
final radio variable sample are denoted as red "+". Objects with complex morphology are 
denoted as blue "x". Objects that are unresolved or marginally resolved by the FIRST are 
denoted as green "x". Objects that fulfill the definition criteria but have been excluded 
from our sample due to contamination by nearby component (s) or unrelated sources are 
denoted as "+" with a box envelope in corresponding color box (see the text for detailed 
description). Note that, due to confusion effects in the NVSS, the asymmetric distribution 
of these quasars around VR = 1 and this asymmetry is more serious for lobe-dominated 
sources than unresolved and marginally resolved sources. The flux density variation VR is 
plotted against radio- loudness RL for our sample of 151 radio variable quasars in the right 
panel. Normal quasars are denoted "+" and BAL quasars by solid circles. 
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Fig. 2. — The SDSS composite quasar spectrum fits (red curve) to the spectra of eight radio 
variable BAL quasars (or candidates labelled with question masks). 
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Fig. 3. — Left panels: the SDSS composite quasar spectrum fits (red curve) for four BAL 
quasar spectra, whose radio emission is resolved by FIRST. Right panels: FIRST survey 1.4 
GHz images of the corresponding BAL quasars; contours are spaced by factors of 2 from 0.5 
mJy. 



